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Abstract. We consider the incidence of Hi absorption in intrinsically small sub-galactic sized extragalactic sources 
selected from sources classified as Gigahertz Peaked Spectrum (GPS) and Compact Steep Spectrum (CSS) sources. 
We find that the smaller sources (< 0.5 kpc) have larger Hi column densities than the larger sources (> 0.5 kpc). 
Both a spherical and an axi-symmetric gas distribution, with a radial power law density profile, can be used to 
explain this anti-correlation between projected linear size and Hi column density. Since most detections occur 
in objects classified as galaxies, we argue that if the unified schemes apply to GPS/CSSs a disk distribution for 
the Hi is more likely. The most favoured explanation for the compact sizes of the GPS/CSSs is that they are 
young sources evolving in a power law density medium. For the GPSs with measured expansion velocities, our 
derived densities are within an order of magnitude of those estimated from ram-pressure confinement of the lobes 
assuming equipartition. Our results therefore support the youth model. 
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1. Introduction 

Gigahertz Peaked Spectrum (GPS) and Compact Steep 
Spectrum (CSS) radio sources have their radio emission 
confined to a very small region (< 10 kpc). Their radio 
spectra and radio morphology indicate this class of sources 
are intrinsically small in contrast to the relativistically 
beamed core-jet sources. 

In these two classes of sources the turn-over frequency 
is anti-correlated with the source size. This has led to the 
argument that that the turnover could be due to syn- 
chrotron self-absorption (Fanti et al. 119901 O'Dea & Baum 
nW7\ Snellen et al. pOOO). A turn-over in the 100 MHz 
regime is characteristic of the CSSs with sizes 1-10 kpc, 
while if the peak flux density instead is seen close to 
1 GHz, the sources are smaller (< 1 kpc) and are clas- 
sified as GPSs. The most likely explanation for the small 
sizes of these objects is that they are young sources which 
will eventually evolve into larger scale classical doubles 
like Cygnus A (e.g. Fanti et al. |1995 Rcadhead et al. 
119961 Owsianik & Conwav il998|l . They therefore provide 
interesting laboratories for studying AGN evolution, and 
present data is consistent with the idea that the sources 
expand with constant velocities while their radio power 
declines as where LS is the linear size (e.g. Fanti 
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et al. I1995|l . The drop of the radio power as the source 
grows in size could indicate that the GPS sources may 
evolve into the lowest luminosity FRII radio galaxies or 
perhaps into the FRIs (O'Dea & Baum 1997). 

The probable youth of GPS/CSSs implies the possibil- 
ity of studying the birth of radio sources; little is known 
of the mechanisms triggering the radio activity. Mergers 
could provide a way to transport gas to the centre of the 
host galaxies and thus be involved in the onset of the nu- 
clear engine. Indeed, optical observations have shown that 
many of the GPS/CSS sources are in disturbed or inter- 
acting systems (de Vries et al. PtiOO O'Dea et al. I1996|l . 
in which we may expect dense nuclear environments. 

A number of studies of gas in GPS/CSS sources have 
already been made at different wavebands. For instance, 
optical studies show strong highly excited line emission 
with large equivalent widths (Gelderman & Whittle fl994fl . 
consistent with interactions between the radio source and 
the interstellar medium. Radio polarisation observations 
have shown that some GPS/CSS sources display large 
rotation measures exceeding > 1000 radm^^ which in 
turn imply large foreground electron densities and mag- 
netic fields (Kato et al. [TWl Stanghellini et al. [TM^ . 
However, the majority of small GPS sources are weakly 
polarised, making it hard to estimate rotation measures 
(e.g. StangheUini et al. I1998|l . The low polarisation sug- 
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gests strong depolarisation, again consistent with a very 
large central density. X-ray absorption observations of 
a few GPS quasars indicate an absorbing column of a 
few 10^^ cm~^ associated with the quasars (Elvis et al. 
I1994|l . However, there exist no extensive study of this phe- 
nomenon in the GPS/CSS class as a whole. Other evi- 
dence for a dense environment comes from free-free ab- 
sorption observations, for example in OQ208 (Kameno et 
al. 2000). Marr et al. I|2()()l|l found evidence for free- free 
absorption in the GPS source 0108-1-388, consistent with 
a 100 pc radius disk with an electron density of 500 cm^''. 
In 1946+708 multi-frequency continuum studies also show 
indications of free-free absorption concentrated towards 
the core and inner parts of the counter-jet, again suggest- 
ing a disk or torus origin (Peck et al. 119991. Disk-like dis- 
tributions of gas have also been found in the optical, and 
the best example so far is the HST dust disk observed 
in the GPS source 4C 31.04 (Perlman et al. EOOIll- Due 
to the sparse number of observations, a typical molecular 
gas mass of GPS/CSSs is not yet known, but CO observa- 
tions of 4C 12.50 have shown a large total gas mass around 
lOi°M0 (Evans et al. lT^ Mirabel et al. lT^g^ . However 
this source is unusual in many other of its properties and 
might not be representative of the GPS/CSS class as a 
whole. 

Another way to study the gas content in these sources 
is by spectral absorption experiments. Such absorption ob- 
servations have the advantage over emission that the sen- 
sitivity is independent of the distance, depending only on 
the background flux density. In absorption experiments 
detection depends on having a significant column density, 
rather than a large gas mass, so this method is sensitive 
to dense gas on small scales in AGN. Using the 21 cm line 
of atomic hydrogen in absorption it is possible to study 
the neutral gas content of GPS/CSS sources. The Hi is 
likely to be only a fraction of the total gas present, how- 
ever such observations provide lower limits to the total 
gas mass and density. The strength of GPS/CSS sources 
at cm wavelengths makes them good targets for such ex- 
periments. In addition their small sizes indicate that lines 
of sight to these sources will sample the dense gas confined 
within the centre of the host galaxy. Similarly the line of 
sight to a GPS source will trace gas within the narrow-line 
region (NLR). 

Here we discuss the results of Hi absorption studies 
in such sources, compiled from our own observations pre- 
sented in Vermeulen et al. H2()03|l . and from the literature. 

2. The sample and its properties 

We are interested in studying the Hi gas content of the 
AGN host galaxy; both the general interstellar matter 
(ISM) in addition to circumnuclear gas disks if present. 
In Tabled we have compiled a list of sources searched for 
Hi absorption which belong to the GPS/CSS class and 
are therefore likely compact with total intrinsic sizes < 10 
kpc. 



In order to achieve a sample of small sources searched 
for Hi absorption which is as complete and unbiased as 
possible, we have selected sources from the following ma- 
jor GPS/CSS lists: Spencer et al. IfWm , the F anti CSS 
sample (Fanti et al. nM1|l . O'Dea & Baum (fTM7| . the CSS 
classified objects from Morganti et al. (1997), de Vries et 
al. <|1997a|l and the COINS sample (Peck et al. lMKTji . The 
selection criteria were: 

1) The redshift must be 0.05 < z < 0.75. This ensures 
a high possibility of the objects having been searched for 
Hi absorption, given the frequency coverage of the WSRT 
(which has a UHF system covering the Hi 21 cm line up to 
z — 0.75). 2) The declination should be larger than —30°; 
the majority of GPS/CSSs are found at redshift between 
0.1 and 1, restricting observations to the frequency flexi- 
ble Northern hemisphere telescope WSRT. 3) The sources 
should have intrinsic linear sizes < 10 kpc, in order to 
represent a true GPS/CSS object. We have therefore ex- 
cluded a few objects (JlllO-f 4817 and J1148-h5924) which 
are known to have a clear one-sided core-jet structure and 
are therefore likely foreshortened by projection (Peck et 
al. I2000|l . After the initial publication from which we se- 
lected the sources, 4 objects (J0432-I-4138, J0702+3757, 
J1227+3635 and J1520+2016) were found to actually have 
redshifts >1.0 and thus have been excluded from our anal- 
ysis. 

This gives a sample of 49 genuine GPS/CSS objects, 
which are listed in TableQ] We find that usable Hi absorp- 
tion data have been taken for 41 of these, most of which 
can be found in Vermeulen et al. H20()3|l . 

Table n]hsts some properties of the sources included in 
our analysis: the source J2000 and B1950 name, a possi- 
ble other name, redshift and the Hi column density with 
the corresponding reference. All column densities assume 
a spin temperature of 100 K and a covering factor of 
unity. Quoted upper limits are given as a 3(T limit, as- 
suming a line width of 100 km s~^. For the sources where 
the absorption line is best fitted by multiple Gaussians, 
the number given is the total integrated column density. 
The FWHM and peak opacity are listed for the compo- 
nent with the largest column density if multiple lines are 
present. The projected linear size was (if possible given the 
source radio morphology) derived from measuring the size 
between the two hotspots in VLBI continuum images. For 
J0556-0241 we have only an upper limits of the size, how- 
ever it is still included in the following analysis. Scant is 
the observed flux density at the redshifted Hi frequency. A 
Friedmann cosmology with Hq = 75 km and qo = 0.5 
is used throughout this paper; all numbers taken from the 
literature are, if necessary, recalculated to this cosmolog- 
ical model. 

3. Results 

3.1. Detection rate 

Including all 41 targets with usable Hi data listed in Table 
n the Hi absorption detection rate is 54%. We show in 
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Table 1. Cols. 1-3; source name, Cols. 4-11; optical redshift, optical host galaxy ID, total observed Hi column density 
or a 3ct upper limit, reference to the Hi absorption data, FWHM and line peak opacity of the detected lines (if multiple 
lines are present, the numbers given are for the component with highest column density), linear size and finally total 
observed flux density at the frequency of the redshifted 21cm line. 



Source 




Other 


2 


Id 


log(ArHi) 


Hiref. 


FWHM 


'^pcak 


Size 


'S'cont 


J2000 


B1950 


name 






cm~^ 




km s~^ 


10^2 


kpc 


Jy 


J0022+0015 


0019-000 


4C 00.02 


0.305 


G 










0.220 




J0025-2602 


0023-263 


OB-238 


0.322 


G 


20.4 


1 


126 


0.93 


2.429 


8.17 


JOlll+3906 


0108-^388 


OC314 


0.66847 


G 


21.9 


2 


94 


44.00 


0.033 


0.18 


J0119+3210 


0116+319 


4C 31.04 


0.060 


G 


21.0 


3 


153 


3.70 


0.091 


2.60 


J0137-H3309 


0134+329 


3C48 


0.367 


Q 


19.3 


4 


100 


0.10 


1.124 


20 


J0141-H1353 


0138+136 


3C49 


0.621 


G 


20.0 


1 


35 


1.39 


5.001 


3.78 


J0224-H2750 


0221+276 


3C67 


0.3102 


G 


<20.1 


1 






9.585 


3.90 


J0348-H3353 


0345+337 


3C93.1 


0.243 


G 


<20.1 


1 






1.186 


2.95 


J0410-f7656 


0403+768 


4C 76.03 


0.5985 


G 


20.4 


1 


61 


1.40 


0.689 


6.30 


J0431-^-2037 


0428+205 


OF 247 


0.219 


G 


20.5 


1 


297 


0.46 


0.653 


4.56 


J0503-h0203 


0500+019 


OG003 


0.58457 


Q 


20.8 


2 


62 


3.60 


0.055 


1.60 


J0542-H4951 


0538+498 


3C 147 


0.545 


Q 


<19.5 


1 






2.717 


28.75 


J0556-0241 


0554-026 




0.235 


G 


<20.8 


1 






<2.420 


0.54 


J0650-H6001 


0646+601 


OH 577.1 


0.455 


Q 










0.014 




J0713-H4349 


0710+439 


01417 


0.518 


Q 










0.118 




J0741-H3112 


0738+313 


01363 


0.635 


Q 


<20.0 


1 






0.041 


1.33 


J0909-H4253 


0906+430 


3C216 


0.670 


Q 


20.1 


1 


177 


0.38 


9.349 


6.19 


J0943-0819 


0941-080 




0.228 


G 


<20.1 


1 






0.148 


3.19 


J1035-I-5628 


1031+567 


OL553 


0.45 


G 


<20.1 


1 






0.163 


1.85 


Jllll-hl955 


1108+201 


OM214 


0.2991 


G 










0.066 




J1120-hl420 


1117+146 


4C 14.41 


0.362 


G 


<19.8 


1 






0.306 


2.75 


J1124-H1919 


1122+195 


3C258 


0.165 


G 










0.242 




J1206-I-6413 


1203+645 


3C 268.3 


0.371 


G 


20.3 


1 


101 


1.00 


5.641 


3.34 


J1252-K5634 


1252+565 


3C 277.1 


0.321 


Q 


<19.8 


1 






6.146 


2.27 


J1308-0950 


1306-095 


OP-10 


0.464 


G 


<20.1 


1 






1.710 


4.82 


J1326-h3154 


1323+321 


4C 32.44 


0.37 


G 


19.9 


1 


229 


0.17 


0.247 


5.41 


J1347-H217 


1345+125 


4C 12.50 


0.122 


G 


20.6 


5 


150 


1.38 


0.166 


4.70 


J1357+4354 


1355+441 




0.64 


G 


21.5 


1 


367 


5.00 


0.067 


0.50 


J1400-K6210 


1358+624 


4C 62.22 


0.431 


G 


20.3 


1 


170 


0.61 


0.218 


5.47 


J1407-H2827 


1404+286 


OQ208 


0.07658 


G 


20.3 


1 


256 


0.39 


0.010 


0.79 


J1414-H4554 


1412+461 




0.186 


G 










0.079 




J1415-^-1320 


1413+135 


OQ122 


0.24671 


Q 


21.0 


6 


18 


34.00 


0.120 


1.25 


J1443-H7707 


1443+773 


3C 303.1 


0.267 


G 


<20.1 


1 






6.295 


2.12 


J1546-H0026 


1543+005 




0.55 


G 


<20.0 


1 






0.049 


1.77 


J1609-H2641 


1607+268 


OS 211 


0.473 


G 










0.223 




J1643-I-1715 


1641+173 


3C 346 


0.161 


G 










8.834 




J1815-H6127 


1815+614 




0.601 


Q 


20.6 


1 


118 


2.03 


0.061 


0.97 


J1816-h3457 


1814+349 




0.2448 


G 


20.7 


7 


80 


3.50 


0.123 


0.80 


J loZO-h / yoo 


loZD-r / yD 




n 90/1 


n 
\j 




1 
1 






u.uoz 


n Q9 


J1829-H4844 


1828+487 


3C380 


0.692 


Q 


<19.3 


1 






6.314 


19.49 


J1944-I-5448 


1943+546 


OV573 


0.263 


G 


20.7 


1 


315 


0.86 


0.149 


1.95 


J1945-H7055 


1946+708 




0.101 


G 


21.5 


8 


357 


5.00 


0.063 


1.01 


J2022-H6136 


2021+614 


OW637 


0.227 


Q 


<19.6 


1 






0.031 


1.94 


J2052+3635 


2050+364 




0.354 


G 


20.9 


1 


16 


16.11 


0.242 


4.51 


J2137-2042 


2135-209 


OX-258 


0.635 


G 


<20.0 


1 






0.914 


4.86 


J2151-H0552 


2149+056 


OX 082 


0.74 


G 


<21.2 


2 






0.016 


0.5 


J2250+1419 


2247+140 


4C 14.82 


0.237 


Q 


<20.1 


1 






0.630 


2.00 


J2344-H8226 


2342+821 




0.735 


Q 


<19.8 


1 






1.006 


5.02 


J2355-H4950 


2352+495 


OZ488 


0.2379 


G 


20.5 


1 


82 


1.72 


0.199 


2.24 



Hi absorption reference: 1) Vermeulen et al. 1200311 . 2) Carilli et al. (119981 . 3) van Gorkom et al. (I1989II . 
4) Hagiwara pers. comm.: tentative detection not yet confirmed, 5) Mirabel ()1989|l . 6) Carilli et al. I|1992^ . 
7) Peck et al. 120001 and 8) estimated from the total integrated VLA spectrum. Peck et al. 119991 
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Table 2. Relation to optical identification. 



10 



22 





Detected 


Nondetected 


Total 


Galaxies 


17 


12 


29 


Quasars 


5 


7 


12 


Total 


22 


19 


41 



Table 3. Relation to radio identification. 





Detected 


Nondetected 


Total 


GPS 


17 


10 


27 


CSS 


5 


9 


14 


Total 


22 


19 


41 



Table 121 that 59% of the galaxies have an Hi detection, as 
opposed to only 42% of the quasars. The latter number 
is however more uncertain since there are fewer quasars 
(29%) in our sample than galaxies (71%). Interestingly, 
we find an even stronger dichotomy when we separate out 
the GPS from the CSS sources by using a linear size divi- 
sion at 1 kpc. Table|21shows that the smaller GPS sources 
have a detection rate of 63% as compared to 36% in the 
larger CSS sources; again the latter number is more un- 
certain due to the smaller number of objects. Below, we 
will discuss at length how this difference could reflect the 
fact that the more compact sources have a larger part of 
their continuum emission covered by nuclear gas. This ef- 
fect has already been suspected on account of the high 
detection rate of Hi absorption in nearby GPS/CSS ob- 
jects fConwav ll996|l . First, rather than simply looking at 
the detection rate in two bins, we show that there may be 
a more continuous relationship between Hi column density 
and linear size. 



3.2. Relation between Nm and linear size 

There appears to be an anti-correlation between the source 
linear size and the largest component of the Hi column 
density as plotted in Fig. (for the objects with multi- 
ple lines we plot the component with the largest column 
density). A visual inspection suggests that the column 
densities for sources < 0.5 kpc are larger than for those 
> 0.5kpc. The data set contains upper limits both in Hi 
column density, as well as a few points which have upper 
limits on their size. We therefore investigated the possible 
correlation in Fig. fusing survival analysis, taking into ac- 
count also the upper limits. We used the software package 
ASURV Rev 1.2 (Lavalley et al. lTM^ . which implements 
the methods for bivariate problems presented by Isobe et 
al. ((1986j . A Kendall's Tau test shows that there exists 



102 



1019 



Detections 
— ^ Upper limits 

< Upper limit size and 

0.1 



0.0' 



Linear size (kpc) 



Fig. 1. Absorbed Hi column density versus projected lin- 
ear size. There is an anti-correlation between the source 
size and the amount of absorbing gas, confirmed by sur- 
vival analysis. Best least square fit taking into account 
the upper limits is shown with the dotted line. This can 
be compared with a least square fit to the detections only, 
plotted with the solid line. 



a correlation between the column density and the linear 
size with a probability > 99%. 

In order to parameterise the correlation we fit the 
simplest model of a power law to the data. The cluster 
of sources with column densities about 10^"'^ cm~^ sug- 
gests that other models are possible, perhaps a two com- 
ponent model or a model with different power laws for 
small and large sources. However, given the uncertain- 
ties of the present dataset we cannot differentiate between 
these models. Linear regression applied with ASURV finds 
the relationship Nm = 7 .2xW^ LS-^-^^cra-'^-, where LS 
is the linear size in kpc. This linear fit is plotted with 
a dotted line in Fig. ^ The solid line represents a least 
square linear fit taking into account the detections only 
(parameterised by Nm = l.^^xW^LS'^-'^^cm-'^). 

While the majority of the upper limits do not show any 
dependence with linear size, there are two of the most com- 
pact sources which have higher upper limits. This raises 
the question whether the observed anti-correlation be- 
tween A'hi and LS could be an observational effect. Since 
the sensitivity of the Hi column density depends on the 
strength of the background continuum, we plot in Fig. El 
the distribution of background continuum strength as a 
function of linear size. Fig.|21shows that we are more sen- 
sitive to small amounts of Hi absorption for the larger 
sources. Clearly, there are no reasons why we should not 
have been able to detect sources in the top right corner 
of Fig. n if they existed, and so we are confident that the 
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— ■ Detections 
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< Upper limit size and 



100 r 



0.1 



_l I 



0.01 0.1 1 

Linear size (kpc) 
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Fig. 2. Observed continuum flux density at observed UHF 
frequency versus projected linear size. This plot shows 
that the lack of large sources with substantial absorption 
is not due to insufficient sensitivity. 



top right hand corner on Fig. ^ is really empty, ft is how- 
ever possible that there are additional absorbers amongst 
the compact sources, with column densities similar to the 
larger sources (~ 10^° - 10^^ cm~^). This would reduce 
the amplitude of the correlation (however the upper limit 
data and its possible dependence on size has already been 
taken account of by ASURV when making its confidence 
estimate) . 

An interesting question is whether the population of 
small sources have larger column densities because the 
FWHM of the absorption line is wider, or because the 
line is deeper. Fig. |3| shows the FWHM versus the pro- 
jected linear size. There is no evidence that there is any 
correlation between the two (tested by a Kendall's Tau 
correlation test). Fig. 0] shows the peak optical depth ver- 
sus linear size. A Kendall's Tau test shows that those two 
variables are correlated with a probability > 99%. The 
solid line represents a least square linear fit to the data 
points. The slope is —0.36, which implies that the slope 
in Fig. n mainly depends on a difference in the observed 
opacity and is not due to differences in the line width. 

3.3. A temporal change in the gas mass or a geometric 
effect? 

There are at least two plausible explanations for the ob- 
served anti-correlation between source linear size and Hi 
column density. The first possibility is that we are tracing 
a temporal change in the total gas mass, and hence in the 
column density. In the evolutionary models of GPS/CSSs 
the smallest sources are the youngest, so we may postulate 




0.1 1 
Linear size (kpc) 

Fig. 3. FWHM versus linear size. There is no correlation 
between the two variables. 



that the younger sources have used less fuel during their 
lifetime, as compared to the older and larger sources. To 
see if this is plausible, it is interesting to consider the to- 
tal gas mass and the rate of use for larger radio sources 
(e.g. FRls and FRHs). The radio loud sources are usually 
associated with elliptical hosts (e.g. Martel et al. 119991 
Bahcall et al. 119971 Dunlop et al. l2()()3|l . which perhaps is 
surprising since it might be expected that more luminous 
sources would reside in more gas rich systems. However, 
the amount of gas needed to sustain an AGN during its 
lifetime is not necessarily large. Consider for instance a 
typical FRI radio galaxy with a bolomctric luminosity of 
the order of 10'*^ erg s^^. Given that these large sources 
have expected lifetimes of ^ 10^ - lO^yrs (Parma et al. 
CnSni), and assuming an energy conversion efficiency of 
10%, we would expect a fuel usage of less than 10^ M0. 
Roughly 15% of the nearby normal ellipticals searched 
for Hi emission display Hi masses > 10* Mq (Huchtmeier 
I1994f) . Amounts smaller than lO"^ - 10* M© are below 
present detection limits, and are thus hard to detect in 
emission. Specifically in ellipticals which harbour power- 
ful AGN molecular gas masses of around 10* - 10^ Mq 
have been found in perhaps 40% of the sources searched 
during CO emission observations (Leon et al. I2()()l|l . with 
detection hmits of the order of 10* M© (Lim et al. 12000(1 . 
We therefore would expect total gas masses of GPS/CSSs 
to be similar to that of classical FRI/FRII sources, given 
the above mentioned expected fuel usage of ^ 10^ Mq and 
that the typical total gas mass of FRI/FRIIs is > 10* Mq. 
Hence, we reject the possibility of the anti-correlation be- 
ing due to a temporal change in the total gas mass, since 
there is no evidence that the radio source would use up 
a significant fraction of its total gas mass during its first 
lO'^ years. 
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1000 r , 

Projected linear size 



100 




Fig. 4. Peak optical depth versus linear size. A clear cor- 
relation is found with a probability > 99%, and the slope 
can explain most of the correlation in Fig. ^ 



Instead we prefer a second explanation, in which the 
observed anti-correlation is a geometric effect of probing 
densities which decrease with radial distance from the cen- 
tre. Here we assume that the sources (which all are < 10 
kpc) are fully embedded within the host galaxy, and that 
the absorbing medium has basically the same distribution 
in all sources (either a spherical or a disk distribution; 
for a cartoon see Fig. O and Fig. If the radio source 
is centrally located in such a gas distribution, the lobes 
of the smaller sources must be embedded in (or in the 
case of a disk, behind) denser gas. Since most of those ob- 
jects are lobe-dominated we thus expect to probe larger 
foreground column densities in the most compact sources. 
Constraining this density profile is useful for constrain- 
ing both the feeding mechanism and evolution model of 
GPS/CSSs, and with this radial density model in mind 
we will in the following section investigate the Hi absorp- 
tion characteristics. 



4. Distribution of the absorbing gas 

For some time it has been suspected that GPS/CSS 
sources have a high Hi absorption detection rate (Conway 
I1996|l . This seems to be confirmed in our analysis of the 
available data (see SectionlHl) . Most of those detections are 
so far done using instruments of modest resolution, which 
means we cannot easily determine where the Hi absorb- 
ing gas is located with respect to the continuum. In the 
following discussion we will first consider the case of the 
absorbing gas in our sample being in a spherical distribu- 
tion, then we discuss the possibility of the Hi absorption 
instead arising in a disk. 




Line of sight 

Fig. 5. Cartoon of the GPS/CSS geometry used for the 
modelling of a spherical density distribution. The radio 
source is confined within the host galaxy, for which we 
assume a maximum radius of 50kpc, i.e. the ISM diameter 
is much larger than the largest projected linear size of 
lOkpc. 

4.1. Spherical distribution 

The strong radio emission in the GPS/CSS sources is most 
often concentrated in two lobes, while the core emission 
(if present) is weak. A cartoon of a GPS/CSS source ge- 
ometry is shown in Fig. O The two lobes do not extend 
beyond the optical host galaxy, hence the probability of 
tracing host galaxy ISM should be high, especially if the 
central environment is enriched due to a merger event. 
Assuming a spherical, radially declining density distribu- 
tion of the ISM gas, the larger sources will probe lower 
density gas. In order to see if such a density distribution 
could explain the anti-correlation in our data, we calcu- 
late the expected integrated absorption TVhi = / n{r)dl 
towards the two lobes, where dl corresponds to the path 
of absorption along the line of sight. In this estimate we 
use a source with an inclination of 45°, however, for this 
symmetric model the resulting Hi column density is only 
weakly dependent on the inclination. 

A first approach is to assume that the density follows 
a King profile: 



= no(l-t-^)-^ 



(1) 



where is the core radius of the gas distribution and 
no is is the density at 1 kpc. From X-ray data, which traces 
the ionised gas component, the core radius for ellipticals 
is usually derived to be '--^ 1 kpc (Trinchieri et al. 119861 
Forman et al. 11985,1 . We consider the simplest case, that 
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Fig. 6. Cartoon of the GPS/CSS geometry used for the 
disk modeUing (not to scale) . The disk model is dependent 
on the disk opening angle 6 and the viewing angle a, see 
Sect.lO 



King's profile 

Radial profile 

Disk model 




10 10 10 10 10 

Linear size (kpc) 

Fig. 7. Results of calculating the expected Hi column den- 
sity given different density profiles. The triangles repre- 
sents the real data. The King density profile, plotted with 
a dotted line, cannot reproduce the same Hi column den- 
sity distribution for any value of /3. Instead, using a power 
law drop of the density provides a more similar distribu- 
tion, with the closest fit for (3= 1.40 and tiq = 9x10^'^ 
cm~'^ (dashed line). Any kind of disk model can fit the 
data reasonably well, here a disk opening angle of 20° is 
plotted (sohd hue). 



the atomic gas component is a constant fraction of the 
ionised gas and follows a similar profile. The two variable 
parameters to be constrained are tlq and (3, and we test 
the goodness of the fit by calculating the sum of squares 
for the observed and calculated Hi column densities. The 
best fit to the data is then achieved for no — 0.05 cm~'^ 
and /3 = 1.06. The resulting absorbed Hi column density 
as a function of source size is plotted with the dotted line 
in Fig. 13 (calculated for a core radius of 1 kpc and host 
ISM outer radius of 50 kpc). In this figure the observed 
values in our sample are plotted with triangles. The King 
model has a more or less uniform density in its centre, 
thus for any value of (3 such a density model is unable 
to reproduce the slope in the data at small linear sizes. 
The only possibility to reproduce the observed data is to 
reduce the core radius by several orders of magnitude and, 
in addition, to increase the core density above 100 cm^^. 
However, then the density model does not follow a typical 
King ISM profile any longer. 

A King model of the density appears inconsistent with 
not only our data, but also with what is generally assumed 
for the medium in which the GPS sources expand. More 
often a simple, steep power law is used, and such power 
laws successfully explain the detection rates of sources 
of different sizes (e.g. Fanti et al. Ji)95; Readhead et al. 
119961 O'Dea & Baum [TIj97|l . Therefore we instead assume 
a power law density profile 

n^n^{L)-P (2) 



where no is the number density (in cm~'^) at ro = 1 
kpc. We again construct a sum of squares surface, 
which gives a best fit to the data for (3 ~ 1.40 and 
no = 0.009 cm^^^tjjjg ^^^^^ jg plotted in Fig. 0). 
Assuming a volume filling factor of unity, the integrated Hi 
mass within a radius of 10 kpc can then be estimated to be 
7x10^ M0. We note that if the filling factor is less than 
1, the estimated mass will decrease, and our estimated 
Hi mass is thus an upper limit. Since the host galax- 
ies of GPS/CSSs may preferentially be found in merg- 
ing/interacting systems (de Vries et al. |2DL)D [ O'Dea et al. 
I199()f) . it is interesting to compare the mass of Hi with 
what has been found in the more nearby population of 
the luminous infra-red galaxies (ULIRGs). Those sources 
are also associated with interactions (although between 
gas rich galaxies), and more or less all of them are de- 
tected in Hi emission with masses in the range of 5x10^ 
to 3x10^" Mq, with typical values of the order of 10^ M© 
(Mirabel & Sanders Il988(l . Even though our objects may 
also be merging systems, a spherical distribution indicates 
that the GPS/CSSs are less gas rich than the ULIRGs. 

Given that the host galaxies of radio loud sources are 
ellipticals, it is also interesting to compare our derived gas 
mass with what is found in nearby non-AGN ellipticals. 
Around 15% of those normal ellipticals are detected in 
Hi emission with masses Mhi ^ 10^ Mq (Huchtmeier et 
al. 119951 Knapp et al. I1985|l . Typical upper limits of the 
remaining 85% are of the order of 10^ - 10^ Mq. Given the 
sensitivity of the Hi emission experiments, the Hi mass we 
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calculate assuming a spherical power law density profile 
is consistent with our GPS/CSS sources being hosted by 
normal ellipticals. 

4.2. Disk distribution 

Instead of spherically distributed gas, another plausible 
distribution is a disk in a plane roughly perpendicular to 
the radio source. For most source orientations this means 
that only the far radio lobe is occulted. It has been shown 
from VLA and ATCA Hi imaging that early-type E and 
SO galaxies often have their Hi in disk systems (e.g. Sadler 
et al. I2UUUI Oosterloo et al. 11999(1 . Many large scale radio 
sources show Hi absorption that is consistent with circum- 
nuclear gas, extended on scales of a few pc in radio loud 
FRIs up to scales of a few 100 pc in low luminosity Seyferts 
(e.g. van Langevelde et al. I2()()()l Taylor 119961 Gallimore 
et alUnni Vermeulen et aL TOI^ . HST imaging has re- 
vealed nuclear disks of gas and dust in radio loud AGN 
(Verdoes Kleijn et al. 199i de Koff et al. lMKHl . In partic- 
ular recent HST observations of a few GPS/CSS sources 
have shown them to contain central disks of gas and dust 
(4C 31.04, 1946-^708 and 1146-^596; Perlman et al. 2001). 
Other high resolution absorption studies of nearby GPSs 
show that the Hi data are consistent with a rotating ring 
(e.g. in 1946-^708; Peck et aLlMiand 4C 31.04; Conway 
I199()|l . Here we consider the case of the absorbing gas in a 
disk oriented perpendicular to the radio axis. 

In this case absorption will only be seen against the far 
lobe. For a cartoon image of the simple disk see Fig.|S| The 
maximum radius of the disk, rdisk, is assumed to be at least 
10 kpc and thus for sources of projected linear sizes up to 
10 kpc the far lobe is always covered by the disk. When 
calculating the Hi column densities, we assume the disk 
has a radial density profile n — no{r /rQ)~^ , where uq is 
the density ai vq — 1 kpc. The calculated column densities 
will be a function of both viewing angle as well as the disk 
opening angle. In our modelling we limit ourselves to a 
viewing angle of 45°, and uses disk opening angles of 10°, 
20° and 30°. We again test the goodness of fit using the 
sum of squares method. Given the scatter in the data all 
such opening angles give similar fits, and in Fig.Qwe plot 
the best fit for a disk full opening angle of 20°, yielding 
no = 0.11 and /3 = 1.45. 

A lOkpc radius disk will, assuming there is an inner 
radius of the disk of Ipc, yield an atomic mass of the 
order of 10^ M0. This is comparable to the results of the 
spherical distribution fScct. HT|l . Following the discussion 
in Sect. 14.11 we find that also for a disk distribution the 
GPS/CSSs appear to be less gas rich than the ULIRGs. 

5. Comparison with GPS/CSS models 

Currently the most favoured explanation to the small sizes 
of most GPS and CSS sources is that they are young, as 
estimated both from their measured expansion speed as 
well as from synchrotron ageing (e.g. Owsianik & Conway 
ITMHl Murgia et al. 11999,1 . In this 'youth' scenario the 



dense gas interacts with the radio plasma giving a large 
radio luminosity which then explains the relatively larger 
number of sources with small size. In order to explain the 
number density of sources of a given size, the radio jet is 
assumed to propagate in a medium with decreasing den- 
sity, approximated with a power law profile oc with 
a slope /3 = 1.5 - 2 (Fanti et al. [TUMI Begelman ITMl 
Readhead et al. ll996|l . Given our uncertainties, it is inter- 
esting that this exponent is comparable with the /3 ~ 1.4 
we get from our Hi absorption. 

Alternatively the GPS/CSSs have been suggested to be 
old 'frustrated' sources prevented from growing by a very 
dense medium (e.g. van Breugel et al. HMl Fanti et al. 
I1990|l . In Sect, n] we mentioned several observations show- 
ing the presence of gas in GPS/CSSs. Other evidence for 
dense environments comes from effects thought to derive 
from interactions with the surrounding medium, several 
CSS objects for example display distorted radio morpholo- 
gies (e.g. van Breugel et al. I1984|) and alignment effects 
between the radio axis and the optical emission (e.g. de 
Vries et al. ll997b|l . Obviously, there is compeUing evidence 
for dense gas in GPS/CSS sources; the question is whether 
the central density is large enough to give frustration. 

While the frustration model is now less popular it is 
still possible that a minority of the GPS/CSS sources are 
created in this way. Here we discuss our Hi results in con- 
text of those two models. 

5.1. The frustration scenario 

In order to confine a median luminosity CSS, an average 
number density in the range 1-10 cm~^ at r = 10 kpc is 
needed (de Young IT993(I . This implies, for a uniform den- 
sity, gas masses of the order of 10^° M©. In our power law 
density model we derive atomic gas masses of the order 
of 10^ Mq, two orders of magnitude less. Comparison of 
molecular data with Hi data for elliptical galaxies have in- 
dicated a Mh2/A/hi ratio of the order of 0.4 (Wiklind et al. 
I1995|l . If such a number is applicable also to the GPS/CSS 
sources, it implies total gas masses around a few times 
10^ M0. This result may change if the ratio A/h2/A^hi is 
very different. For the gas rich ULIRGs however, the to- 
tal molecular gas mass derived from CO emission, and 
the atomic gas mass derived from Hi emission are of the 
same order (10^ M0; Downes & Solomon 119981 Mirabel 
& Sanders Il988|l . and will thus change the resulting total 
gas mass by a small factor. Given those estimates, the Hi 
absorption data suggests that there is not enough gas to 
frustrate the radio sources. 

5.2. The youth model 

The young ages estimated from the hotspot advance 
speeds are strong arguments in favour of the youth model. 
Using ram pressure arguments, the lobe advance speed v^, 
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Table 4. Density of the medium where the lobes expand, as expected from lobe advance speed (next), compared to 
Hi density estimates from a spherical density model (ns.iooK) and a disk model (?t.d,ioook)- Estimates are done using 
h = 0.75. The first part of the table lists sources which are detected in Hi absorption, while the second part lists 
sources that are not detected in Hi absorption. 





Source 


radius 




Pi 






Chi 


WS.IOOK 


wd.ioook 






kpc 






cm"^ 




cm"^ 


cm"^ 




Sources with 


JOlll+3906 


0.017 


0.098 


1.1x10"" 


1 


3.6 


9.1 


24.6 


368.4 


detected Hi absorption 


J1407+2827 


0.005 


0.050 


2.5x10"" 


2,3 


31.8 


0.1 


1.5 


23.9 




J1944+5448 


0.087 


0.210 


2.5x10"" 


4,5 


1.8 


1.1 


0.3 


4.2 




J2355+4950 


0.078 


0.101 


9.7x10"^^ 


6,7 


3.0 


0.7 


0.2 


3.1 


Sources without 


J0713+4349 


0.058 


0.126 


1.9x10"" 


8,7 


3.8 




0.5 


6.8 


detected Hi absorption 


J1035+5628 


0.082 


0.310 


2.9x10"^ 


9 


0.1 




0.3 


4.1 




J2022+6136 


0.011 


0.060 


6.4x10"" 


10 


56.5 




5.0 


76.1 



^References for lobe advance speed and the lobe pressure: 1) Owsianik et al. 119981 2) Polatidis & Conwav 120031 
3) Stanghellini et aLQW! 4) Polatidis et al. il9Ml 5) A. Polatidis, pers. comm., 6) Owsianik et al. 119991 
7) Conway et al HWI 8) Owsianik & Conwav IT^^ 9) Taylor et al. 1^17001 and 10) Tschager et al. 1^17001 



is related to the external mass density poxt and the hotspot 
pressure pi by 



For the spherical density profile, we use the power law 
profile (Eq. |5J)) and express the density ns as a function 
of radius and spin temperature: 



Pi = Pcxtvl = ricxtmpvl 



(3) 



ns 



CHinHi(^)=CHinor,-^,(^) 



(4) 



Thus, for sources with measured expansion velocities, 
we can estimate the density rioxt at the hotspot location. 
We can further compare next with the density expected 
using the power law density distribution models discussed 
in Sect.H 

Published in the literature we find a selection of 7 
sources with measured advance speeds and quoted (or 
sufficient information to calculate) internal hotspot pres- 
sures. For those sources, we list in Table ^ the distance 
between the core and hotspot, the advance speed, hotspot 
pressure and the corresponding references. We use Eq. © 
to derive external densities next (Table 0} at the given 
hotspot radius. 

In order to derive densities using the disk and spheri- 
cal power law density profile (Eq. |(21)), we first note that 
all sources listed in Table 01 have very small linear sizes. 
Especially, for sources < 0.1 kpc there are large offsets be- 
tween the measured Hi column density and the Hi column 
density as predicted by our simple models (Fig.0). For the 
sources with measured Hi column density, we therefore ap- 
ply a correction factor (Chi = ^hi, measured /^hi, predicted), 
in order to better estimate the density in the source in 
question. We also note that the estimated density will de- 
pend on the spin temperature; the values for the column 
density used in Sect. assumes a T^p = lOOK. However, 
close to an X-ray irradiating AGN the spin temperature 
is likely to be above lOOK, and probably closer to 8000 K 
if the gas is purely atomic (Maloney et al. .1996)1 . 



Applying the best fitting values of tiq and (3 from 
Sect. O and using a spin temperature of lOOK we 
achieve: 



ns,iooK ^ Chi 0.009ri7p^, 



(5) 



At a given radius Eq. Q gives the Hi density, and is 
thus a lower limit to the total density. Using the radius 
of the sources in Table 0] we calculate the corresponding 
values of the density ng.iooK- For this model we have cho- 
sen to use Tsp — lOOK, since that gives the closest match 
to the densities ricxt estimated from the ram-pressure ar- 
gument. If the spin temperature is 8000 K, the estimated 
ng will be 80 times higher. 

If instead the absorbing gas is in a disk distribution, 
we can estimate the density at the hotspots (which must 
be in a medium external to the disk, see Fig. by as- 
suming that the disk internal pressure will be similar to 
the pressure outside the disk. Since we do not have any 
information on the properties of such a disk, we will first 
consider the simplest case, assuming the disk to be com- 
pletely atomic and to have an uniform kinetic temperature 
Tdisk — Tsp. The external medium is assumed to be ionised 
with an electron temperature Text = 10* K. Applying the 
ideal gas law the external density no can be written as: 



nD = ^HinHi(^)(^) = CHinor,-/e( 



lOOK^^Text' 



(6) 
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For a disk opening angle of 20°, we derived a value of 
no ~ 0.11 cm^'^ and /3 — 1.45 (Sect. and assuming a 
spin temperature of lOOOK we rewrite Eq. (jSJ into 

riD,ioooK ^ 0.11 CHiri7pi/5 (7) 

Again, this is a lower limit to the total gas density 
present. In order to get densities which best agree with 
those estimated via the ram pressure argument, we here 
use a Tsp = lOOOK. The resulting values of riD.ioooK are 
hsted in Table H 

Excluding JOlll-l-3906, for the spherical power law 
model most next and ng.iooK values are remarkably sim- 
ilar, and in most cases only differ within an order of a 
magnitude; which is consistent with the GPSs propagat- 
ing through a purely atomic spherical gas distribution. 
The offset can depend on a number of things. Firstly, it 
is possible that the spin temperature has a very different 
value. Secondly, and closely linked to the spin tempera- 
ture, is the expected molecular abundance. If the kinetic 
temperature is as low as lOOK, the molecular abundance 
is expected to be high (Maloney et al. I1996|l . and thus 
our Hi density estimate is only a lower limit to the total 
gas density. Thirdly, the value of P and no could be dif- 
ferent, given the sparse amount of data used in the fitting 
(22 points). Within a 90% confidence region the value of 
P may change as much as ±0.25 and no as much as ±4e-3, 
with a corresponding difference of 1 in the number den- 
sity at a radius of 0.05 kpc. Finally the individual viewing 
angles of the sources in Table 0] may affect the assumed 
radius where the hotspot is located. 

The densities derived from the disk profile are equally 
similar, differing with only a factor of a few from n^xt 
(except for JOlll-l-3906). The offsets are not surprising, 
given the large number of uncertainties in the calculation 
of n-D.ioooK- In addition to the reasons already given for 
the spherical distribution, the ratio T^isk/Tcxt could for 
instance be very different from the assumed value. 

We conclude that within the scope of the power law 
profile, we more or less can reproduce the densities ex- 
pected in the youth model for the GPS sources with mea- 
sured lobe velocities. 

6. Orientation effects on the observed HI 
column density: galaxies versus quasars 

In most cases the spatial resolution of present Hi data 
is not enough to determine whether the absorption cov- 
ers both lobes (consistent with ISM gas) or only one lobe 
(indicating a disk distribution). In the previous sections 
we showed that, given the limited information we have 
about the density profiles, both a disk as well as a spher- 
ical distribution are consistent with the observed relation 
between Hi column density and linear size. One possibil- 
ity to distinguish between those models is to look at any 
possible orientation effects. Assuming the unified scheme, 
quasars are supposed to be seen at a smaller viewing angle 
a (more end-on) than the radio galaxies. In fact, we have 
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Fig. 8. Hi column density versus fraction of core emission 
fc- Quasars are marked with crosses, and galaxies with 
squares. In general the quasars are found in sources with 
larger core prominence, consistent with smaller viewing 
angles. 

shown that there appears to be a larger detection rate in 
galaxies than in quasars (Sect.|OJ- 

For a given projected source size of 1 kpc, the amount 
of expected Hi absorption as a function of viewing an- 
gle can be estimated. Given that for a projected size of 
1 kpc we have a Hi column density sensitivity of around 
lO^*' cm^^ (Fig-D; the spherical density profile predicts 
we would detect equal column densities in galaxies (with 
say a ^ 45°; Barthcl 1989) and quasars (a ^ 45°). In 
contrast for the disk profile there will be a range of opti- 
mum viewing angles, for which we are likely to detect Hi 
absorption. This exact range of viewing angles will be ex- 
tremely dependent on the disk model used. In Sect. 16. ll we 
discuss if the galaxies and quasars in our sample are con- 
sistent with having different viewing angles, as expected 
from unified schemes. If so, a spherically symmetric distri- 
bution could be considered less likely, since objects classi- 
fied as galaxies appear to have higher Hi column densities 
than quasars. 

6.1. Core prominence 

Because of relativistic boosting the core strength is consid- 
ered to be a good indicator of orientation; the more dom- 
inant the core is, the smaller the viewing angle. For larger 
radio sources, this has been investigated by a number of 
people (e.g. Orr & Browne HW^ . Since the GPS/CSS 
objects in general have weak or no core emission, simi- 
lar studies have been difficult. Stanghellini et al. H2001|l 
reports on VLBA observations of 10 GPS sources. The 
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authors conclude that in the few cases where cores were 
detected, the core of the GPS galaxies accounted for ~ 2% 
of the total flux density, while in the GPS quasars the core 
accounted for ~ 20%. Moreover, Saikia et al. (jl995|il stud- 
ied a sample of CSS sources with detected radio cores, and 
found that the degree of core prominence were consistent 
with those for larger radio sources and also consistent with 
the unified scheme. In the study by Saikia et al. (jIP95|, 
there are 7 sources that can also be found in our sample 
(3 quasars and 4 galaxies). For the galaxies the average 
fraction of emission from the core is — 0.012, while 
the corresponding number for the quasars is = 0.26. In 
Fig. m we plot the Hi column density versus /c for those 
objects, and we mark quasars and galaxies with different 
symbols. Clearly the quasars have larger core prominence. 
This implies the quasars have a smaller viewing angle. 

Only one of these quasars is detected in Hi absorp- 
tion, while three of the galaxies are Hi absorbers. For the 
sources in common with the sources of Saikia et al. (4995)) , 
it appears as if the upper limits of the Hi column density 
for the quasars are smaller than the detections for the ra- 
dio galaxies. It would be interesting to carry out deeper 
observations for a larger sample of the CSS quasars with 
upper limits; if deeper searches showed (for a majority of 
the sources) lower column densities or upper limits for the 
quasars than for the radio galaxy detections, this would 
strongly argue for a disk origin for the Hi absorption. 

6.2. GPS quasars versus GPS galaxies 

It should be noted that it is debated whether the GPS 
quasars are really intrinsically small sources like the 
GPS galaxies. Snellen H1997|l proposes that GPS galax- 
ies viewed end-on will have a flat spectra and thus the 
GPS galaxies and quasars are different types of objects. 
This has some support in a very different redshift distri- 
bution for the GPS galaxies and quasars (Snellen 1997ii 
Stanghellini 1199^ . and the fact that many GPS quasars 
appear to have core-jet morphologies when observed at 
high angular resolution f Stanghellini .2001|) . If this is the 
case, the low detection rate in the GPS quasars reported 
here could be explained by them being large scale radio 
sources, hence with a lower probability of the line of sight 
going through the host galaxy. 

7. Implications for the spectral turn-over 
mechanism 

The integrated radio spectrum turn-over frequency for 
GPS/CSS sources is known to have a relationship i^t oc 
^5^-0.65 (o'Dea & BaumTWTl). Most often this is ascribed 
to synchrotron selfabsorption (e.g. O'Dea & Baum 1997) : 
however an alternative possibility is the one of free-free 
absorption (e.g. Bicknell et al. 1997). A few VLBI results 
indicate the presence of a nuclear disk causing free-free 
absorption (Marr et al. lMITl Peck et al. ITU^ . 

Since our Hi results show that the column density 
along the line of sight depends on the source linear size 



(and thus turnover frequency), it is likely that a similar 
relationship exists for an ionised gas component. The free- 
free absorption turn-over frequency is related to the elec- 
tron density, vt oc (n^)^/^'^. If we simply assume the gas 
fraction to be constant {nm/ric = C) and A^hi — n^L 
where L is the absorbing pathlength, the observed rela- 
tionship (A^Hi oc LS~^''^^) implies a free-free absorption 
caused turnover relationship Vt LS'^° '*^. This is much 
less steep than observed {h't oc LS*""-^^), implying that 
free-free absorption can not be the only turn-over mecha- 
nism present. 

Conversely we could assume that the observed rela- 
tionship between turnover frequency and linear size is 
due to free free absorption. Then, vt oc LS~'^'^^ implies 
that the slope for the ionised gas component necessarily 
is larger than the slope we observe for Hi. In this case 
perhaps there is a gradient of the fraction of ionised gas 
decreasing with radius. 



8. Conclusions 

We have investigated the combined results of Hi absorp- 
tion in GPS/CSS sources. The most striking result is that 
the smaller GPS sources (< 0.5 kpc) tend to have larger Hi 
column densities than the larger CSS sources (> 0.5 kpc). 
This anti-correlation between linear size and absorbed Hi 
column density is mostly an effect from opacity differences 
rather than line width differences. It is possible to use a 
simple power law profile of the density to explain the ob- 
served Hi column densities, and both a spherically sym- 
metric medium as well as disk can be fitted to the data. 

Our modelling cannot distinguish between the spher- 
ical and disk density models, and most observations at 
present do not have enough resolution to determine the 
location of the absorbing gas. However, to date high reso- 
lution Hi data exist for a limited number of sources that in 
all cases indicate Hi disks fConwav ll99fil Peck et al. 119991 
Peck & Tavlor [r998|l . This argues for a disk model, and 
in addition we notice that most Hi absorbers are found in 
sources classified as galaxies. By studying our sources in 
their radio continuum properties we find that the galax- 
ies and quasars in our sample behave similarly to larger 
scale sources, where orientation effects are known to take 
place. Therefore, it is easier to argue for a disk model, 
where sources with very small viewing angles in addition 
to sources observed edge-on would be less likely to have 
Hi detectable in absorption against their lobes. The fact 
that quasars have smaller column densities argues against 
a spherically symmetric model for the Hi gas. 
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